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FINAL  REPORT 

USARO  Contract  No.  DAAG29-80-C-0001 
January  1,  1980  through  February  28,  1983 

Shock-Tube  Kinetics  Studies  of  Formaldehyde 
Oxidation  Reactions  Involving  N02»  NO,  and 
N2O  With  and  Without  O2 

K.  G.  P.  Sulzmann 
Principal  Investigator 
Energy  Center 

University  of  California,  San  Diego 
La  Jolla,  CA  92093 

I .  Statement  of  Problem 

The  detailed  combustion  chemistry  of  solid  propellants  in  use  by  the 
U.S.  Army  is  not  understood.  Of  particular  interest  are  the  chemis¬ 
try  of  nitrocellulose  as  well  as  double  base  and  nitramine  propel¬ 
lants  in  general,  and  the  gas-phase  reactions  of  cyclotetramethy-  j 

lene-trinitramine  (RDX)  and  of  cyclotetramethylene  tetranitramine  1 

1—3  I 

(HMX)  in  particular.  It  appears  that  reactions  between  aldehyde  I 

and  nitrogen  dioxide,  besides  nitric  oxide  and  nitrous  oxide,  domi-  I 

nate  the  gas-phase  combustion  processes  near  the  propellant  surface  j 

in  these  systems,  and  it  is  important  to  obtain  a  quantitative  1 

understanding  of  the  overall  kinetics  for  engineering  applications  1 

and,  in  the  long  run,  a  detailed  quantitative  knowledge  of  the  domi-  I 

nant  reaction  mechanisms  and  their  relevant  rate  coefficient  for  1 

future  computor  modeling  of  the  gas-phase  propellant  chemistry.  I 

Some  shock-tube  measurements  of  importance  for  these  systems  I 

1-3  I 

have  been  conducted  in  the  past.  However,  they  were  based  upon  I 

in  situ  measurements  of  only  one  fuel-intermediate  (CHgONOg  or  1 

CHgCHgNOg  in  ref.  3}  or  only  an  oxidizer  species  (NOg  In  refs.  1  I 
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and  2  or  NgO  in  ref.  2),  and  product  analysis  was  performed  after 
quenching  of  the  reacting  mixture.  Since  shock- tube  quenching 
methods  usually  cannot  yield  reliable  information  about  the  degree 
of  reaction  progress  at  a  certain  time  after  initiation,  additional 
in  situ  shock-tube  measurements  were  required,  particularly  for 
the  fuel  species  besides  intermediates  and  products.  Specifically, 
the  oxidation  of  formaldehyde  by  NOg,  besides  NO  or  NgO,  are  of 
importance  in  this  respect. 

Previous  studies^  at  relatively  low  temperature  (390  <  T/°K  <  465) 
and  at  about  0.4  atm  total  pressure  showed  that  the  reaction  in  homo¬ 
geneous  and  first  order  in  both,  the  formaldehyde  and  the  nitrogen 
dioxide  concentration  with  an  overall  rate  coefficient  whose  pre-expo¬ 
nential  factor  and  whose  apparent  activation  energy  increased  with 
increasing  temperature.  These  studies  were  conducted  in  a  static 
reactor  and  the  pressure- increase  was  used  as  a  measure  of  the  reac¬ 
tion  progress  with  the  result 


-(H2C0)’=  -0.724(N02)’  *  k(H2C0)(N02)  , 


where  K  -  1.26  x  1010  exp(-  15,000/RT)  cm3/mole-sec  for  391  <  T/°K  <  433 

and  E  *  1012  exp(-  29, 000/RT)  cm3/mole-sec  for  433  <  T/°K  <  465. 

1  2 

In  subsequent  shock-tube  studies  ’  in  argon-diluted  mixtures  of 
HgCO  and  N02  at  elevated  temperatures  and  pressures,  the  reaction 
progress  was  monitored  by  spectral  absorption  of  N02  in  the  visible 
and  resulted  in  an  overall  rate  law  for  NOg-consumption  of  the 
approximate  form1 


-CN02)*  * 


(2) 


CO) (M)  (N02) 

with  Ep  *  1.26  x  io13  exp(-  26,700/RT)  cm3/mole-sec  for  992  <  T/°K 
<  1828  and  7.3  <  p/atm  <  12.3.  Thus,  the  dependence  on  reactants 
as  well  as  the  overall  rate-coefficient  appeared  to  change  with 
increasing  temperature  and  pressure. 

Consequently,  we  have  conducted  additional  measurements  behind 
reflected  shock-waves  for  temperatures  of  835  <  T/°K  <  1175  and  for 
total  pressures  of  0.96  <  p/atm  <  2.85  for  varying  initial  concen¬ 
trations  of  H2CO  and  N02  by  using  HeXe-laser  absorption  at  3.508u 
for  monitoring  R2CO  as  well  as  infrared  emission  near  3.5y  and  near 
7 . 5y  for  H2CO  and  N02- 

In  the  following  we  summarize  important  results  which  were  ob¬ 
tained  during  the  last  two  years  of  this  research  effort,  after  re¬ 
direction  of  our  studies  on  the  chemistry  of  combustion  inhibition  by 
CF^Br  after  the  first  contract  year. 

II .  Summary  of  Most  Important  Results 

A.  Overall  Rate  Parameters 

We  have  obtained  empirically  overall  induction  (I),  fast  re¬ 
action  (R),  and  equilibrium  (E)  times  for  our  experimental  condi¬ 
tions.  The  induction  time  tj  was  defined  by  the  time  interval  be¬ 
tween  the  reflected-shock  arrival  at  the  laser  observation-station 
and  the  intercept  of  the  maximum  slope  for  the  HeXe-laser  trans¬ 
mission  rate  with  the  initial  transmission  I  after  reflected  shock- 

o 
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heating.  Our  data  could  be  fitted  best  by  the  relation 


ti(h2c)o(N02)o/(M)o 


(H2C0)o(N02)o 

(H2CO)- 


(3) 


where 

£j  =  3.03 xio17  (±76%)  exp[-29,100(+5.2)/RT]  cm3/mole-sec  (3a) 

is  an  effective  rate-parameter,  and  the  quality  of  the  least-squares- 

2  - 
fit  is  r  *  0.941.  The  initiation  rate-parameter  kj  is  approximately 

related  to  the  rate-coefficient  k^  for  the  initiation  step  of  the 

H2C0-N02-reaction  by 

k  -  rt  -  Un  W  ,3b) 

1  ■  I  '  (In  I„/I0)  (M)0  •  (3b) 

where  Ia  *  laser  transmission  after  reaction  completion,  ID  *  initial 

—2 

laser  transmission,  and  Ij  *  laser  transmission  at  t * Tj,  where  F *  10 
(see  below).  The  experimentally  established  relation  (3)  shows  that 
the  induction-rate  is  first-order  in  both,  the  initial  H2CO-  and  the 
initial  N02-concentration. 

The  equilibrium  times  rE  where  defined  by  the  time  interval  be¬ 
tween  the  ref lected-shock  arrival  at  the  laser  observation-station 
and  the  intercept  of  the  maximum  slope  for  the  HeXe-laser-transmission 
rate  with  the  final  transmission  las  after  completion  of  the  reaction. 


•  J  i“ 
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The  data  could  be  fitted  best  by  the  relation 


-  ,  <H2C0)o(N02>o 

te(H2C0)0(N02)o/(M)o  I 


where 


kE  *  2.03  * 1016(±85%)  exp [-  22,800  (±7.6%) /RT]  cm3 /mole-sec  (4a) 


is  the  effective  rate-parameter  for  reaching  equilibrium  with  a 
least-squares-f it  quality  of  r3  =  0.892.  The  empirical  result  (4) 
shows  that  the  overall  rate  for  equilibration  is  also  first  order 
in  both  the  initial  HgCO-  and  the  initial  NOg-concentration.  Thus, 
it  can  be  implied  that  the  overall  oxidation  of  HgCO  with  N02  is 
first-order  throughout  in  both  (HgCO)  and  (N02). 

From  the  fit-data  for  tj  and  te  followed  for  the  fast-reaction 


tin>e  tR  te"TR) 


(h2G0)  (n°2)0 

VH2C0>o<N02V<M>o  =-  -  TH2C0)— 


where 


£„  «  4.02  x  1015(±13%)  exp [-  18,100  (±1.6%)/RT]  cm3 /mole-sec  (5a) 


is  the  effective  rate-parameter  for  the  rapid-reaction  regime  at 

2 

constant  pressure  and  the  least-squares-f it  quality  is  r  ■  0.998. 
The  rapid-rate  parameter  ER  is  approximately  related  to  the  overall 
rate-coefficient  kR  for  the  maximum  HgCO  consumption- rate  by 


(5b) 


kR  =  ERGr 


(H2C0)o/0I)o  r 
a.+  io)  an  ijiQ) "  kR 


with  the  same  meaning  of  symbols  as  given  under  equation  (3b)  and 
-1  -2 

Gr  *  10  -  10  (see  below). 

From  the  results  (3)  through  (5)  follows  that  the  consumption 
of  H2CO  by  oxidation  with  N02  is  first-order  in  HgCO  and  first-order 
in  N02  as  was  also  found  in  reference  4,  but  not  by  the  observations 
in  references  1  through  3.  The  apparent  activation  energy  decreases 
from  29.1  kcal/mole  during  the  induction  regime  to  18.1  kcal/mole 
during  the  fast  oxidation  regime.  Here  the  activation  energy  for 
the  rapid  oxidation  regime  is  similar  to  that  (-19  kcal/mole)  ob¬ 
served  above  433°K  in  reference  4,  although  the  pre-exponential  for 
the  effective  rate-coefficient  kR  is  larger  than  that  for  E  by  fac¬ 
tors  of  10-100. 


B.  Initiation  Step 

Rate  data  for  the  initiation  step  were  obtained  by  measuring 
the  early  HgCO  consumption  rate  at  the  time  of  reaction-initiation 
immediately  after  shock-heating.  The  HeXe-laser  transmission  ex¬ 
hibited  a  finite  initial  rate  in  all  cases.  Analysis  of  the  rates 
in  terms  of  the  initial  rates  for  HgCO  resulted  in 


-(H2C0) 


ki(H2C0)(N02)  . 


(6) 


Here, 


kl  -  (M)o*(N02)0 


-i,.„  .-i  r(1,v*-> 


r<in  i 

L(1»  iL/r.)tJt<0 


3.46  X  1016  (±102%)  exp [-29 , 600( ±7 . 1% ) /HTJ  cm3 /mole-sec  (6a) 


is  the  rate-coefficient  for  the  initiation-step  from  a  ? east-squares- 

fit  with  a  quality  of  r2  *  0.970. 

The  result  demonstrated  that  the  initiation-step  is  first-order 

in  both  the  H2CO-  and  the  NOg-concentration,  as  expected  from  the 

—2 

relation  (3)  above,  and  it  confirms  that  F  =  10  in  equation  (3b). 
Within  the  experimental  errors,  the  activation  energy  E^  for  the 
initiation-step  (6)  is  the  same  as  the  effective  activation  energy 
Ej  for  the  induction-regime  (3a). 

The  relatively  large  pre-exponential  factor  and  low  activation 
energy  for  k.  can  be  rationalized  by  invoking  unimolecular  rate 

5 

theory  in  the  high  pressure  limit  for  the  formation  of  a  vibration- 
ally  excited  complex  Q*  =  (N02  •  HgCO)  which  is  subsequently  decom¬ 
posed  by  collisions  with  a  third  body  (M)  to  products  R^  and  R2 
according  to 


H2C0  +  N02  t  Q*  , 

(i) 

Q*  +  M  -*•  Rx  +  R2  +  M  . 

(d) 

Here,  some  or  all  of  the  products  can  react  further  with  H2C0  and 


NOg  by  chain  and/or  regular  reactions.  For  (Q*)  in  a  steady  state 

and  k^M)  »  follows  the  initiation-step  (6).  The  rate-coeffi- 

5 

cient  kj  can  be  obtained  from  a  "primitive"  unimolecular  theory, 
viz. , 

k  -  io13e-Ei/RT  <Ei/RT)f~1  (6b) 

Ki  “  ( f-1 ) ! 

which  satisfies  (6a)  with  a  reasonable  number  of  vibrational  degrees 

*  -1 
of  f  *  5-6  excited  in  Q  and  a  reasonable  decay  time  of  t_^  =  k_i  » 

—13  * 

10  sec  for  Q  .  It  is  interesting  to  note  that  this  expression 
extrapolates  within  the  experimental  uncertainties  to  the  values  of 
E  obtained  for  the  slow  initial  reaction  at  lower  temperatures  in 
reference  4.  Thus,  it  appears  that  the  same  initiation  step  applies 
at  low  and  at  elevated  temperatures. 


C.  Fast  Oxidation  Regime 

Rate  data  for  the  rapid  oxidation  regime  were  obtained  by 
determining  kR  directly  from  the  relation  [comp.  (5b)] 


<1.  “ 


V 


o  tR 


(* 


GrER) 


(7) 


which  applies  for  the  maximum  rate  of  Hg CO-consumption  with  Im  the 
HeXe-laser  transmission  at  the  time  of  maximum  HgCO-consumption.  A 
least-squares-f it  to  the  data  results  in 


-  7.85  x  1013(±122%)  exp[-  16,000  (+14.7%)/RT]  cm3 /mole-sec  (7b) 


J 


R 


10 
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with  a  quality  of  r  *  0.780.  The  activation  energy  for  kR  is,  with¬ 
in  the  experimental  uncertainties,  the  same  as  that  determined  for 

-1  2 

kR  in  equation  (5a),  and  the  expression  (7b)  confirms  Gr  »  10  -  10 

14  3 

in  the  relation  (5b).  The  pre-exponential  factor  of  -10  cm  /mole- 
sec  in  (7b)  is  reasonable  for  a  binary  reaction  rate-coefficient. 

It  should  be  noted  that,  within  the  error  limits,  ER  =  E^/2  which 
implies  a  short  chain  mechanism  for  the  rapid  reaction  regime  with 
chain  branching  and  chain  propagation  steps  of  low  or  negligible  acti 
vation  energies  and  with  a  bimolecular  chain-termination  step  involv¬ 
ing  two  chain  carriers.  Since  kR  was  determined  with  T  taken  as  the 
temperature  at  the  initiation  of  the  reaction,  and  since  T  increases 
during  the  course  of  the  oxidation  process  by  potentially  highly 
exothermic  reactions,  the  heat-release  will  influence  the  maximum 
oxidation-rate  and  the  value  for  ER  above  is  a  lower  limit  for  the 
actual  overall  rate  coefficient  for  the  maximum  oxidation  rate. 
However,  this  effect  is  small  in  view  of  the  dilution  of  the  system 
with  total  reactant  concentration  of  $15%,  and  amounts  to  less 
($+12%)  than  the  experimental  error  in  ER. 

Elucidation  of  the  detailed  reaction  mechanism  requires  further 
study,  particularly  in  view  of  the  differences  in  the  reactant  orders 
for  the  overall  oxidation  process  obtained  by  us  as  opposed  to  those 
references  1  through  3  and  because  of  our  confirmation  of  its  first 
order. 


D.  Additional  Results 

We  have  obtained  evidence  for  the  existence  of  a  formaldehyde 
dimer  in  the  gas-phase  at  room  temperatures  and  at  high  temperatures 
after  shock-heating  in  pure  and  argon-diluted  mixtures  of  formalde- 

g 

hyde.  The  line-strengths  for  the  monomer  (m)  and  dimer  (d)  of 

HgCO  near  the  HeXe-laser  wavelength  of  XL  *  3.508u  as  a  function  of 

-2  -1 

temperature  are  (in  cm  -atm  )  for  285  <  T/°K  <  1600,  respectively, 

vT> =  s£<vT>ra  (8) 

with  m  =  1.94  (i7.3%)  and 


Sd(T)  =  (SK/b)°(T/T0)n+*exp[+(AH0/RT0)(T0/T-l)]  (9) 


with  n  =  1 . 07 (-37%) ,  (SK/b)°  *  3.67  ( ±0.  7%)  cm-1- atm-1,  and  AHQ  »  940  cal/ 

7 

mole  a  theoretical  estimate  for  the  heat  of  formation  for  the  dimer  species. 

HeXe-laser  absorption  by  argon-diluted  NOj  at  3.508u  at  room- 

temperatures  and  elevated  temperatures  revealed  no  dependence  on 

total  pressure  and  can  be  expressed  for  285  <  T/°K  K  2000  by  an 

—1  —1  8 

absorption  coefficient  (in  cm  -atm  )  of  the  form 


P(T)  =  A(T„/T)m+i{exp[-(E0/RT0)(T0/T-l)]} 


(10) 


x  [l-exp(-hcwL/kT)] 


where  A  *  1.37  x  io-^  cm-1-atm-1,  m  *  2.54,  and  EQ  -  5.66  kcal/mole 


2 

with  a  least-squares-f it  quality  of  r  =  0.885. 

Approximate  analytical  relations  for  the  spectral  intensity 

distribution  and  the  radiative  recombination  rate-coefficient  for 

g 

Br-atom  free-bound  transitions  have  been  derived  for  Br-atom  con¬ 
centration  measurements  in  the  context  of  initial  studies  on  the 
chemistry  of  combustion  inhibition  by  CF^Br  before  re-direction  of 
the  imphasis  of  our  research  towards  H2CO-oxidation.  In  addition, 
a  simple  method  for  calculating  Voigt-prof iles  by  a  programmable 
hand-calculator  (HP-34C)  was  developed3-0  which  is  useful  for  laser- 
absorption  studies  on  gases. 
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